Shapes of Semifiexible Polymer Rings 
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The shape of semifiexible polymer rings is studied over their whole range of flexibility. Inves- 
tigating the joint distribution of asphericity and nature of asphericity as well as their respective 
averages we find two distinct shape regimes depending on the flexibility of the polymer. For small 
perimeter to persistence length the fluctuating rings exhibit only planar, elliptical configurations. At 
higher flexibilities three-dimensional, crumpled structures arise. Analytic calculations confirm the 
qualitative behavior of the averaged shape parameters and the elliptical shape in the stiff regime. 
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It is a well-known fact dating back to 1934 that the 
shape of a flexible coil is overall prolate From the 
isotropy of space, the intuitive expectation would be a 
spherically symmetric conformation. However, this idea 
implies rotational averaging and in fact entropy is maxi- 
mized for a single trajectory of a polymer if the number 
of segments in each direction is inhomogeneous. 

After a series of theoretical investigations, based on 
both analyses @, H, H| and simulations [B|, Q, only 
with the onset of single molecule techniques, experiments 
could prove the asymmetric shape of a flexible polymer 
0, 0] and address the relevance of a polymer's shape in 
biology. The overall shape of a polymer is important 
for its mobility in heterogeneous media such as cyto- 
plasm and the depletion forces between larger complexes 
in polymer solution [To| . For the transcription of viral 
genome or plasmids, the shape of its DNA might enhance 
or reduce the accessibility for enzymes depending on the 
spatial distance between DNA-segments DNA as 



many biopolymers is semifiexible, behaving like a ther- 
mally fluctuating elastic rod on length scales of the order 
of its persistence length. Considering the shape of viral 
DNA and plasmids this limit is applicable and deserves 
investigation. In fact, most of the short genomes as well 
as plasmids are circular, yielding an even stronger con- 
straint for the polymer's shape. To obtain a complete 
picture for any circular DNA, it is desirable to under- 
stand the shape of semifiexible rings as their flexibility 
is varied. A polymer's shape is well characterized by the 
asphericity [2j as the deviation from spherical symmetry. 
The degree of prolateness or oblateness is captured by 
the independent nature of asphericity || . Measurements 
of both their mean values give a good indication of how 
the average outline looks like, but fail to reflect the total 
ensemble of configurations which can only be accessed 
via the shape parameter's distribution. 

We employ Monte Carlo simulations to study the shape 
of semifiexible polymer rings over a large range of flexi- 
bility To give a complete picture of the polymer's change 
of shape as its flexibility increases, the joint distribution 
of asphericity and nature of asphericity as well as their 



respective averages will be presented. We find two dif- 
ferent shape regimes. In the first, the flexibilities are 
small resulting in dominantly planar polymer ring con- 
figurations. In the second, at large flexibilities, crumpled 
three-dimensional (3D) structures prevail. In both the 
stiff and the flexible limit analytic calculations explain 
the observed behavior. 

Characterizations of the shape of a polymer's trajec- 
tory {r(s)}, s € [0, L], are based on the radius of gyra- 
tion, primarily a measure for the spatial extent. Gener- 
alizing to a radius of gyration tensor Q, 

Qi] = j j J dsri(s)rj(s)-r^ J dsr^s) Jdsrj(s), (1) 

the eigenvalues A, of the tensor describe the spatial ex- 
tent along each principal axis. Measuring the variance 
of the eigenvalues, the deviation from a fully symmetric 
object is obtained, denoted asphericity, A. Furthermore, 
prolateness or oblateness of the object is specified by the 
nature of asphericity, E, measuring the skewness of the 
eigenvalues. Choosing the normalization such that the 
quantities are independent of the total length, the as- 
phericity of a polymer is defined by Q 
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where Qij = Qij — <5yTr Q/3. The nature of asphericity 
is given by [5[ 
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The asphericity takes values < A < 1, where A = cor- 
responds to a spherically symmetric object. For A — 1, 
the polymer is fully extended, forming a rigid rod. The 
nature of asphericity is bounded between — 1 < E < 1. 
E = — 1 is obtained for a fully oblate object such as a 
disc, while E = 1 is the result for a prolate object as 
a rigid rod. As the asphericity and the nature of as- 
phericity are independent, a joint distribution yields a 
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thorough classification of stochastic objects such as ther- 
mally fluctuating polymers. For reasons of comparison, 
we will adopt the parameters p = 2yA £ [0, 2] and 
8 = arccosS/3 € [0,7r/3] for the joint distribution de- 
fined by Cannon et al. [5(. 

Those parameters are directly connected to the eigen- 
values of the radius of gyration tensor by: Ai = A(l + 
pcos(6>)), A 2 = A(l + pcos(8 - 2tt/3)) and_ A 3 = A(l + 
pcos(9 + 2tt/3)), where Ai > A 2 > A 3 and A denotes the 
mean eigenvalue. Using these relations, the shape dia- 
gram presented in Fig. [T] is constructed. In the region 
of both large p and large 8 one eigenvalue becomes neg- 
ative excluding these parameter sets for real structures. 
Along the solid line separating the excluded conforma- 
tions from possible ones, at least one eigenvalue is zero. 
Hence, the solid line represents all planar configurations 
ranging from the fully oblate geometry of a rigid ring with 
6 = 7r/3 via elliptical shapes to the fully prolate structure 
of a rigid rod at 8 = 0. Below the solid line, 3D conforma- 
tions are exhibited as all eigenvalues are now greater than 
zero. The shape is rather oblate for 8 > 7r/6 or compara- 
tively prolate for 8 < n/6 as illustrated by the ellipsoids 
enclosing a polymer's trajectory. Towards smaller p, the 
structure becomes less and less aspherical resulting in a 
spherically symmetric conformation for p = 0. 

For flexible open polymers the shape distribution is 
known to be almost exclusively prolate and highly as- 
pherical being peaked around 8 = ir/40 and p — 1.55 [5] 
as indicated by the diamond in Fig. [T] However, the con- 
formation of a rigid ring lies just at the opposite end of 
the shape diagram at 8 = ir/3 and p = 1. As the states of 
a highly flexible ring polymer can be assumed to be sim- 




9 = arccos E/3 

FIG. 1: (color online). Illustration of the overall shape of 
polymer configurations depending on the asphericity A and 
the nature of asphericity E. Along the solid black line the 
structures are planar; configurations beyond that line are ex- 
cluded as they do not correspond to real structures. The di- 
amond indicates the peak in the distribution of flexible open 
polymers. The sequences of circles and triangles denote the 
mean shape of open and ring polymers at integer flexibilities 
starting from L/l p = at (0, 2) and at (n/3, 1), respectively. 



ilar to those of flexible open polymers, a strong crossover 
between a stiff and a flexible regime seems inevitable. 
Heuristically we may argue that the shape of the tight 
fluctuating ring in the stiff limit is expected to be dom- 
inated by the first modes since higher modes are almost 
not thermally excited. Both the first in-plane "breath- 
ing" mode and the first transverse bending mode yield an 
elliptical conformation as can be illustrated by deforming 
a strip of paper connected to form a ring. Although the 
ring rotates in space when fluctuating, the elliptical shape 
itself remains planar, being oblate for small eccentricities 
and becoming prolate for large eccentricities of the el- 
lipse. Towards the flexible limit also higher modes are ex- 
cited resulting in a crossover to the flexible regime where 
the conformations are three-dimensional and crumpled 
as expected for a closed random walk. 

The METROPOLIS Monte Carlo method was em- 
ployed to simulate a discretized semiflexible ring of to- 
tal length L and persistence length l p . The ring is de- 
scribed as a polygon composed of N tethers of fixed 
length a = (L/ir) sin(7r/7V) and direction t. The energy 
assigned to an individual configuration is given by the 
elastic energy, E = NkBT(l p / £)Xh=i(1 — t,t i+ i), im- 
posing periodic boundary conditions, tjv+i = ti. New 
conformations are achieved by pivot moves fl2j | , perform- 
ing 10 6 Monte Carlo steps per segment. Measured ex- 
pectation values of the mean square diameter (D 2 ) were 
in accordance with analytical expressions (]~3j | up to the 
estimated statistical error. The effect of self- avoidance 
is neglected, as its impact on the shape of even flexible 
polymers was shown to be only of the order of 1% 0]. 

The change of shape as the flexibility increases is best 
studied when analyzing the shape distribution at differ- 
ent flexibilities as plotted in Fig. [2l We distinguish be- 
tween a stiff regime exemplified by L/l p = 1,4 and a 
flexible regime represented by L/l p = 16,32. The ge- 
ometry of a ring induces an apparent stiffening of the 
fluctuating polymer to approximately five times its un- 
constrained flexibility [13(. Therefore, even L/l p = 4 
belongs to the stiff limit. By comparison with the shape 
diagram in Fig. [T] polymer ring configurations in the stiff 
regime are identified to be almost exclusively planar rang- 
ing from totally oblate to comparatively prolate shapes. 
In the flexible regime, crumpled structures that fill 3D 
space dominate the broader configuration space. The dis- 
tribution changes from rimlike being strongly peaked in 
the asphericity to lenslike with the major weight on pro- 
late and highly aspherical conformations, although less 
rodlike than observed for open polymers. In agreement 
with experimental observations [7|, l8j the distribution of 
shapes is very broad yielding also very extended confor- 
mation close to p — 2. 

In between the two asymptotic shape regimes a 
crossover is observed represented by L/l p = 8 in Fig. [2] 
During this crossover both crumpled, 3D configurations 
and planar structures are almost equally probable yield- 
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FIG. 2: (color online). The distribution of asphericity A and 
nature of asphericity E at different levels of flexibility. For 
tight rings, L/l p < 5, planar conformations dominate, while 
the configurations become truly 3D beyond L/l p of the order 
of ten. The crosses indicate mean and variance of each sin- 
gle shape parameter deviating from the most probable state. 
Note the change of the color scaling as the distribution spreads 
out. 



ing the largest spread of well-occupied conformations in 
the configuration space. Beyond the stiff regime in Fig. [5] 
the mean of each single shape parameter deviates from 
the states with largest joint probability, showing its lim- 
itation in identifying a polymer's configurations. Identi- 
fication of the conformations and insight into the width 
and the form of the density of states are only attainable 
by the distribution of the shape parameters, highlighting 



their importance for studying polymer shapes. 

Investigating the mean asphericity and the mean na- 
ture of asphericity, qualitative arguments can be quanti- 
fied and compared to previous results in the limit of infi- 
nite flexibility. The change of both shape parameters on 
increasing flexibility L/l p is depicted in Fig. [3] For a rigid 
ring, the asphericity is given by A — 0.25, being fully 
oblate: £ = — 1. Up to L/l p sw 5 both asphericity and na- 
ture of asphericity grow linearly with the flexibility, obey- 
ing (A) stiff = 0.25 + 0.01 L/l p , (£> stiff = -1 + 0.3 L/l p . 
This linear dependence classifying the stiff regime is ex- 
plained by the shape of an ellipse whose axes grow and 
shrink with the square root of the flexibility, respectively, 
as will be discussed in the next paragraph. A similar scal- 
ing argument has been given by Camacho et al. 1J] ana- 
lyzing planar rings. Beyond this stiff regime a maximum 
of the mean asphericity is reached. Increasing the flex- 
ibility further, higher modes become accessible. These 
undulations contract particularly the major axis of the 
"ellipse" , decreasing the variance of the eigenvalues of the 
radius of gyration tensor and hence yielding a declining 
asphericity. The asphericity approaches the exact value 
for an infinitely flexible polymer ring, a closed Gaussian 
chain, (A)a iC = 0.2464, derived by Diehl and Eisenriegler 
3] in a power law with exponent v — —1.3. Compared 
with a flexible open polymer with (A) flj0 = 0.396(5) [| 
a polymer ring is much more spherical. Analytic calcu- 
lations based on a perturbation expansion of a closed 
Gaussian chain for finite flexibility forecast a positive 
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FIG. 3: (color online). Monte Carlo simulation data for the 
mean asphericity (A) and the mean nature of asphericity (E) 
versus increasing flexibility L/l p . Both grow linearly with 
raising flexibility for tight rings (solid line). Then, (E) sat- 
urates at a prolate shape, while the asphericity decreases in 
a power law (dashed line). Error bars are of the size of the 
symbols. 
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correction in first order [l5(, and, hence, explain why 
the asphericity approaches its Gaussian limit from above. 
Also for two-dimensional polymer rings our analytic ar- 
guments predict a linear increase of the asphericity in the 
stiff limit and a monotonic decrease in the flexible limit 
explaining the nonmonotonic behavior of the shape pa- 
rameter observed in previous simulations in two dimen- 
sions [3, [3] ■ Experiments found a slight increase of the 
asphericity versus L/l p in the flexible regime for linear 
DNA 7| not to be justified with our predictions of the 
shape of semiflexible polymers discarding for the sake of 
simplicity additional DNA effects, such as twist, nicks or 
supercoils. Although DNA rings make up a huge field of 
biological processes where shape matters, resolving their 
3D shapes is a challenge due to their small size and rate 
of change. Ring polymers of larger per sistence length 
l p such as cytoskeletal filaments as in [17( or mesoscopic 
polymer materials are feasible to measure our results. 
The nature of asphericity increases monotonically satu- 
rating at a value of (£)fl, c = 0.53 at L/l p « 15, being less 
prolate than a flexible open polymers with (E)fl iQ = 0.745 
[5| . Towards the Gaussian limit the sequences of averaged 
shape parameters of ring and open polymer as depicted 
in Fig. [T] can neither cross nor depart from each other. 
Therefore, both their (E) approach their limiting value 
monotonically. Overall, the geometric constraint induces 
a bias towards more spherical and oblate structures. 

The linear growth of (A) and (E) for tight rings is ana- 
lytically predictable based on the assumption of a planar 
shape. Due to the first bending modes, the ring of ra- 
dius R c becomes an ellipse, where the major and minor 
axis are the radius R c increased and decreased, respec- 
tively, by yj the amplitude of the undulations of a 
weakly bending rod. In the weakly bending limit fluctu- 
ations parallel to the average axis of the contour are sec- 
ond order to undulations perpendicular, resulting in the 
approximate bending energy E — ^ksTl p J Q ds r" (s) 2 , 

-1/2 



d3/2;- 



l£\m,\lM, where 7 



and yielding yj (r j_ 
denotes a numerical constant. Hence, the asphericity and 
the nature of asphericity are equated in the limit of small 
flexibilities L/l p , where the first modes truly dominate 
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0.25 + 2 1 L/l p + O((L/l p ) 2 ) , (4) 
-1 + 54 7 L/Z P + 0((L/Z P ) 2 ) . (5) 



These analytic results forecast the observed behavior. 

Concluding we have employed the joint distribution of 
asphericity and nature of asphericity and their respective 
means as well as analytic arguments to show that the 
shape of semiflexible polymer rings exhibits two distinct 
regimes depending on their flexibility. Tight rings are 
planar "ellipses", while flexible rings are 3D, crumpled 
structures. These two regimes may have implications for 
a variety of biological processes such as the flow behav- 
ior or the accessibility of DNA rings to enzymes. As the 



shape of stiff, elliptical rings may not be changed by hy- 
drodynamic forces considerably since they will behave as 
rigid discs, rings in the flexible regime may undergo tum- 
bling motion with alternating collapse and stretching as 
observed for flexible open polymers [20j. Similarly, the 
time it takes an enzyme to find its assigned binding site 
on a DNA strand should be larger if the DNA confor- 
mation is planar, as the enzyme cannot easily travel to 
DNA-segments separated afar along the backbone by 3D 
diffusion as in coiled up 3D structures. In this context 
of opposed behavior in the two shape regimes, polymers 
in the crossover region where both shapes are equally 
probable may show striking properties. Depending on 
the manner a polymer changes between a planar and a 
crumpled shape, e.g. randomly or following a particular 
trajectory, and its timescale, a broad variety of biologi- 
cal functionality can emerge. The characterization of a 
semiflexible polymer by its shape can therefore enable a 
coarse-grained modelling of complex biological processes. 

Financial support of the German Excellence Initiative 
via the program "Nanosystems Initiative Munich (NIM)" 
and of the Deutsche Forschungsgemeinschaft through 
SFB 486 is gratefully acknowledged. 



[1 
[2 

[3; 

[4. 
[5 

[6 

[T 

[8 

[9 
[10 

[11 

[12 
[13 
[14 

[IS 
[16 

[17 

[is; 

[19 
[20 



W. Kuhn, Kolloid-Z. 68, 2 (1934). 

J. A. Aronovitz and D. R. Nelson, J. Physique 47, 1445 
(1986). 

H. W. Diehl and E. Eisenriegler, J.Phys. A: Math. Gen. 
22, L87 (1989). 

K. Sole, J. Chem. Phys. 55, 335 (1971). 

J. W. Cannon, J. A. Aronovitz, and P. Goldbart, J. Phys. 

I 1, 629 (1991). 

M. Bishop and J. H. R. Clarke, J. Chem. Phys. 90, 6647 
(1989). 

C. Haber, S. A. Ruiz, and D. Wirtz, Proc. Natl. Acad. 
Sci. U.S.A., 97, 10792 (2000). 

B. Maier and J. O. Radler, Macromolecules 34, 5723 
(2001). 

R. Verma et al, Phys. Rev. Lett. 81, 4004 (1998). 

M. Triantafillou and R. D. Kamien, Phys. Rev. E 59, 

5621 (1999). 

T. Hu, A. Y. Grosberg, and B. I. Shklovskii, Biophys. J. 
90, 2731 (2006). 

K. V. Klenin et al, J. Mol. Biol. 217, 413 (1991). 
K. Alim and E. Frey, Eur. Phys. J. E 24, 185 (2007). 

C. J. Camacho, M. E. Fisher, R. R. P. Singh, J. Chem. 
Phys. 94, 5693 (1991). 

K. Alim and E. Frey, unpublished. 

R. E. Norman, G. C. Barker, and T. J. Sluckin, J. Phys. 

II 2, 1363 (1992). 

M. M. A. E. Claessens et al, Nat. Mat. 5, 748 (2006). 
J. Shimada and H. Yamakawa, Biopolymers 27, 657 
(1988). 

T. Odijk, J. Chem. Phys. 105, 1270 (1996). 

C. M. Schroeder, R. E. Teixeira, E. S. G. Shaqfeh, and 

S. Chu, Phys. Rev. Lett. 95, 018301 (2005). 



